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Refraction at the sea surface 

Comparison between spectral distribution of downward 

( — ) and upward ( ) irradiance for a solar elevation 

of 55-60° in the Sargasso Sea. (Lundgren and 
Hojerslev, Kobenhavns Universitet, Rep. Inst. Fys . 
Oceanog. L4 (1971) 

Compilation of irradiance data by G.L. Clarke and 

G.D. Ewing, Chapter 17 in Ref. 2. Spectral of the 

energy of upwelling light on 13 August 1969 in 

Buzzards Bay at the altitudes shown (measured by 

TRW spectrometer) and at 0*20 m below the sea 

surface (measured by Scripps spectrometer) as 

compared with downwelling irradiance of the NASA 

Standard spectrum outside the atmosphere, and the 

downwelling irradiance measured on 14 August 1969 

at the sea surface. For the "above surface" curve 

see text. The energy scale on the right is for 

downwelling light; that on the left is in relative 

units, except that the values on the left also 
-2 2 

represent 10 W/m /nm for the upwelling light at 
0*20 m. 
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Upwelling spectral radiance just below the surface 
computed from measurements at depths of 2.2 and 7.7 m. 


ii 


Fig. 

Fig. 

Fig. 

Fig. 

Fig. 
Fig. 
Fig . 

Fig. 

Fig. 

Fig. 

Fig. 


7. Component parts of the upwelling radiance signal, L . 

z 

8. Computed apparent spectral radiance of the ocean (and 
its components) as observed above the atmosphere. 
Blue-green water. 

9. Computed apparent spectral radiance of the 
(and its components) as observed above the 
Blue water. 

10. Computed apparent spectral radiance of the 
(and its components) as observed above the 
Green water. 

11. Diffuse reflectance. 

12. Single scatter event. 

13. Apparatus for measuring diffuse reflectance of 
pollutants . 

14. Attenuation curves in the near ultraviolet and in 

the visible part of the spectrum, s is total scattering 
coefficient for pure water and puie sea water as a 
function of wavelength. Data compiled by A. Morel, 
Chapter 1 in Ref. 2. 

15. Attenuation curve for water between 0.2 and -2.8jjm, 
compiled by A. Morel, Chapter 1 in Ref. 2. 

16. Geometry of polarized sky reflection. 

17. Reflectance of undisturbed air/water boundary for 
polarized light (r and r ) and for unpolarized 
light (r) . 


ocean 

atmosphere . 
ocean 

atmosphere . 


Fig. 18. Reflectance and transmittance of water/air boundary 

as a function of the angle of observation for several 
windspeeds. Curves at left are for radiance coming 
up from below the surface; those on the right are 
for radiance incident from above [Ref. 4.]. 

Fig. 19. Sky luminance distributions in sun-zenith plane for 

several locations and solar zenith angles, 0 . Data 
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compiled by Austin [Ref. 3.]. 

Fig. 20. Round-trip attenuation. 


Fig. 21. Vertical beam transmittance of the atmospheric models 
(curves A, B, C, D) and as measured by Guttman 
(curves a and b) . Applied Optics 7 (1968) p. 2377. 


Fig. 22. Absorption cell for measuring pollutant transmittance 
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Fig. 23. Quantities for remote sensing equation, low altitude 
case. 


Fig. 24. Polluted layer over clean seawater. 

Fig. 25. Light from polluted area diluted by scattering. 


THEORY OF SPECTRAL RADIANCE OF POLLUTANTS AT SEA 


This report concerns remote measurement of pollutants dumped 
in the sea, not oil slicks, but soluble pollutants that change 
the color of the water. The sensor is a spectral radiometer 
that flies over the polluted area and compares its spectral 
radiance (color) to that of surrounding clean seawater.. It 
is fairly easy to make a qualitative determination that a 
plume of pollutant occupies a certain area. However, the 
goal here is quantitative, namely to infer the concentration 
of pollutants using the measured radiance of the sea compared 
to laboratory measurements of reflection and transmission 
•spectra of the pollutants. Many previous papers and reports 
have discussed remote sensing of ocean color. The main 
contribution here is how to uses data quantitatively that 
compare the polluted area to surrounding clean water. 

We treat the subject in three steps. The first section defines 
the quantities involved and describes means for measuring 
them. The second section derives the equations for remote 
sensing with a low-flying aircraft, in which case the absorption 
and radiance of intervening air is negligible. The third 
section applies to high-flying aircraft and satellites, in 
which case the radiance of intervening air is the major problem. 
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1.0 


FUNDAMENTALS 


This section describes measurements of light intensity, reflection, 
and transmission. 

1.1 Irradiance 

The brightness of the sea surface is proportional to the amount 

of incident sunlight measured as an irradiance, so much power 

2 

(watts) per unit area (m ) . This incident light is a combination 
of collimated direct sunlight and diffuse sky light. The angles 
of incidence are not too important as long as the sun is fairly 
high in the sky because refraction at the sea surface bends the 
rays toward the nadir, Fig. 1, so that all the light is directed 
more or less downward. 

A meter that measures irradiance (irradiometer) is nothing more 
than a lensless light detector with a precise area A of detection 
defined by a mask. Fig. 2. The detector measures power P, and 
irradiance is given by 

H = P/A (1) 

Since A is constant, a practical meter is calibrated directly in 

2 : 

irradiance units, such as watts/m . A suitable antiref lection 
coating on the detector helps to make the meter insensitive to 
the _igle of incidence. To measure light incident on the sea, 
it is important that the irradiometer be exactly horizontal. 
Obviously it will intercept more sunlight if tilted toward the 
sun. 

Since we are interested in color of the sea, we need a spectral 

2 

irradiometer which measures spectral irradiance, watts/m /nm. 

The irradiometer is converted to a spectral irradiometer by 
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Fia. 2. Refraction at the sea surface 





covering it with a series of absorbing filters of bandwidth B. 
(Interference filters would be sensitive to angle.) Then the 
spectral irradiance is given by 

H x « H/B - P/AB (2) 

Numerous measurements of irradiance have been made and compiled 
in texts on ocean optics, e.g. Jerlov [Ref. 1] pp 71, 128 (Figs. 
63 and 65b) , and 144 and Jerlov and Nielsen [Ref. 2] pp 88 and 
399. Two summaries of irradiance appear in Figs. 3 and 4. Note 
that just beneath the suriace, downwelling ifradiance in the 
blue-green is about 20 times the upwelling. This ratio increases 
to 35 in the red. In a later estimate we use the value 

H(510 nm) = 1.2 W/ra 2 /nm 
1.2 Radiance 

The remote sensor that looks at the sea cannot be a lensless 
irradiometer. An imaging sensor is required to distinguish 
a plume of polluted water from surrounding clean water. The 
amount of light power this sensor detects from each pixel is 
proportional to the solid angle of rays that it intercepts, 
and so the appropriate measure of sea brightness is power 
(watts) per unit area (m ) per unit solid angle (steradian) , 
a quantity called radiance, L. 

Figure 5 shows the optics in- a radiometer. Power received at 
the detector is given by 

/area of resolution cell\ x /solid angle intercepted 
p “ L x \Fig. 5, lower right / \Fig. 5, lower left 

From the geometry, these quantities are evaluated as follows: 


4 


m nm 



). Comparison between spectral distribution 

of downward ( — ) and upward ( ) irradiance 

for a solar elevation of 51,-60° in the 
Sargasso Sea. (Lundgren and Hojerslev, 
Kobenhavns Universitet, Rep. Inst. Fys. 
Oceanog. 14 (1971)) 





eompUatiion of irradiance data by G.L. Clarke 
and G C. Ewing, Chapter 17 in Ref. 2 Spectra 
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where A » HD2/4 is the collecting aperture , and D its diameter, 

f the focal length/ and h the altitude of the sensor. Note 

2 

that h drops out, and Aa/f is a constant of the optics. This 

means that the detector output can be calibrated directly in 

2 

terms of radiance, watts/m. /nm. Also note that the radiance 
of the image is the same as that of the object, neglecting 
absorption in the lens: 


L 


i 



P 

(A/f 2 )a 



L 


(4) 


Again we are interested in color of the s>ia, and so spectral 
radiance is required: 

L x * L/B (5) 

Absorbing filters couli ‘e used, but in this case the optics 
can be used to collimate light in a suitable way to use a 
grating or prism monochromator. These details are beyond the 
scope of this report. From this point on, we drop the subscript 
X with the understanding that all radiances and irradiances 
are spectral. 


Austin [Ref. 3] has compiled graphs of radiance that are 
particularly useful for remote sensing. Four of these appear 
in Fig. 6 and 8 through 10. Fig. 7 defines the notation he 
uses. The following features in these date influence the 
techniques discussed later: 


^ L (510 nm) * 0.9 * 10~ w/n; /sr/nm regardless of 
the color of the water (Fig. 6) 

^ In the deep red, sky radiance reflected by the 
surface (L r ) exceeds that which has penetrated 
the water iL w ) . The crossover falls in the range 
580 - 650 nm. 
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Fig. 6. Upwelling spectral radiance just below 
the surface computed from measurements 
at depths of 2.2 and 7.7 m. 



fig. 7. Component parts of the upwelling 
radiance signal. L . 
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^ From high altitude or outer space, the observed 

radiance (L z ) is dominated by the veiling luminance 
of the intervening air (L*) . 


1.3 Diffuse Reflectance 


Consider direct sunlight incident on a rough surface as shown 
in Fig. 11. The ratio of reflected to incident power is 
reflectance p, a dimensionless quantity. In order to relate 
L seen by the radiometer to H, the incident light, we need 
the ratio. L to H, a quantity with dimensions of per steradian: 


L 

H 



( 6 ) 


where (3^^ is the effective solid angle over which the light 
is scattered. We call this ratio the diffuse reflectance.* 


Consider the bulk reflectance R of a liquid, excluding reflection 
at the surface, which we treat separately. This liquid can 


* In the case of a perfectly white diffuse reflector, p = 1 

gle fo: 

= 1/tt. This relationship 


and ft.** = it, the effective solid angle for a lambertian 


Hence R. 


eff 

(cosine law) surface. 

provides a means for calibrating a radiometer against an- 

irradiometer given a near-perfect surface. For airborne 

instruments, a cloud bank serves this purpose very well 

providing i.t is thick enough to stop direct sunlight (cannot 

see the sun from below) but not thick enough to be an absorber 

(dark thunderhead) . In this case 

L + L , = H/tt (7) 

up down 

where values of L are measured from above the clouds looking 
down, and below looking up anc; for a cloud of the right thickness 
they are about equal. 
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be either clean or polluted seawater; in either case it contains 
particulate matter that reflects light diffusely. It is 
possible in principle to derive an expression for R using the 
attenuation a of the medium and the diffuse reflectance adz, 

Fig. 12, caused by particulate matter (plus Rayleigh scattering) 
in each infinitesimal layer. A simple expression useful later 
applies to the case in which the absorption rate, a, dominates 
over scattering: 

a > s = o(0)du or. a > 4Tra(180°) 

Multiple scattering is negligible in this case, and the single 
scatter rays in Fig. 12 lead to the expression 

R = J Q e (adz) e = a J Q e dz 

R = a/ 2a, a > 4ira(180°) (8) 

It is possible to treat more general cases, but we see no point * 
in relating R to the bulk optical properties of the liquid. 
Instead it is more reliable simply to measure reflectance in 
the laboratory for various concentrations of the polluted 
layer. This is especially true in view of the fact that some 
pollutants undergo chemical change as the concentration changes. 
For example, the acid waste that DuPont dumps precipitates iron 
when diluted. 

1.3.1 Reflectance of Pollutants 


Figure 13 shows an apparatus suitable for measuring diffuse 
reflectance. A known irradiance H illuminates a cylindrical 
column of liquid (thickness z) . A radiometer in a side arm 
measures L, and then 

RT = L/H (9) 
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Fig. 12. Single scatter 
event 


Fig. 13. Apparatus for measuring 
diffuse reflectance of 
pollutants 




Here T is a transmission factor for the air/liquid dielectric 
interface (discussed later) . It is just a well to measure the 
product RT instead of each separately, because the product is 
what appears later in the equations for remote sensing. 

Several features of this apparatus (Fig. 13) are noteworthy. 
First the cylindrical wall below the surface should be highly 
reflecting, silvered preferably on the inside (and coated to 
prevent corrosion) . This simulates the effect of an unbounded 
ocean by replacing light that would come from outside the 
cylinder (dashed ray lower right) by. a reflected ray on the 
inside. Other surfaces should be black to reduce stray light. 
To avoid edge effects, the incident beam must overfill the 
cylinder, and the radiometer's field of view must not extend 
to the side walls. 


1.3.2 Diffuse Spectral Reflectance (Color) of Clean Seawater 

We tend to think that the color of seawater is obvious, somewhere 
between blue and green, but actually the apparent color is 
contaminated with a lot of sky reflection, somewhere between blue 
and white, and the correction is not always straightforward. 
Probably the only reliable means to measure the true color of 
seawater alone, i.e. R w , its spectrum of diffuse reflectance, 
in all weather requires a boat trip to the site in question. 

One measures radiance L u (upwelling) with a radiometer just 
under the surface looking down, and irradiance (downwelling) 
with an irradiometer just under the surface looking up. Then 


R w * 


V H d 


In a few sites where pollutants are dumped regularly, it may be 
practical to do this enough times to establish the variations 
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of R w with season, hour of the day, and sea state. Then, 
from that time on, R w may be considered a known spectrum 
in the remote sensing equation. 

Obviously we need a fail-back position for areas that do not 
justify all that labor. One alternative is to fly over the 
areas and take measurements when the sea is calm and the sky 
uniform. Under these circumstances, the remote sensing 
equation (Eq. 17, discussed later) may be solved for R w since 
all other quantities are- either known or easily measurable. 
One has to hope that the values of R found in this way also 

w 

hold when the weather changes. 


The final fall-back position is 'to use only two colors where " ' 
R^ is fairly invariant from place to place and therefore 
known. These values are 


-3 -1 

Green: R,,(510 nm) = 8 x 10 sr 

w 

Near infrared: R w (700 to 1100 nm) = 0 


( 10 ) 


These numbers require some explanation. The green color is 
where blue water, blue-green water, and green water all have 
the same reflectance. Fig. 6. This assumes that the water in 
question is moderately deep and its color is not influenced 
by bottom reflection, silt, or gelbstoff. Unfortunately we 
know of no case in which R^ (510 nm) or any other wavelength 
has been measured directly. We calculated the value above 
from a ratio L/H in which L and H were typical values measured 
in widely different parts of the world by different groups 
who presumably do not compare instruments. We took L from 
Fig. 6 and H is an average from Figs. 3, 4, and others like 
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them. This situation should be corrected because any quantity 
that is fairly invariant should be pinned down and published 
as a universal semi-constant whose range of variation is known. 

The second of Eqs. 10, R w (700 to 1100) ■ 0, simply means 
that R^ is negligible. To estimate the actual value, let us 
use Eq. 8 with 

-3 

c(160°) = 10. per meter per steradian 
a ■ a » 0.5 per meter. 

The former avoids a fairly sharp peak at 180° and represents 
a fairly large value of the scattering function in a backward 
direction. The value of the attenuation constant a comes 
from Fig. 14. Now Eq. 8 gives 

R ■ o/2a = 10~^ per steradian, 

an order of magnitude less than the green value, which already 
represents a very dark scattering medium. (Deep seawater 
looks almost black compared to most diffuse reflectors.) The 
useful range of infrared cuts off at about 1100 nm where light 
penetrates only about a centimeter of water. Fig. 15, hardly 
enough to see an appreciable amount of diluted pollutant. 

In view of the difficulties in determining R , it is advisable 

w 

to formulate the remote sensing equations in such a way that 
this quantity drops out. This is done in Sec. 2, but the 
result is not usable in all circumstances. 

1.4 Sky Reflection 

Reflections from fog, foam, spray and particulates in the sea 
are all of the diffuse type described in the preceding section. 
There is one exception, reflection at a calm sea surface. 
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This is a simple specular reflection with no angular spread, 
and so collimated direct sunlight is still collimated after 
reflection, and the radiometer sees a reflected image of the 
sun. The brilliance of this image is overpowering and masks 
the color of the sea. Any radiometer looking down at the sea 
must be pointed away from the sun's image to obtain useful 
data. With this precaution, the downlooking radiometer sees 
the color of the sea, but it is still contaminated by a 
significant amount of reflected sky radiance: 


where S is the radiance of the sky as measured looking up at 

the angle supplementary to the direction of L gJc as Fig. 16 

shows. The extent of this contamination is considerable: 

L is about 1/3 of the desired radiance L from the water 
r w 

(Figs. 8, 9, and 10) throughout most of the spectrum, but 

the two cross over in red and the contamination exceeds the 

signal. Since no angular spreading occurs at reflection, 

2 

and both L g ^ and S have spectral radiance units (W/m /sr/nm) , 
it follows that r is a dimensionless number like the reflectivity 
of any mirror surface. 


Values of r are plotted in Fig. 17, which shows that it is 
dependent on polarization, the symbol "j." referring to polarization 
(E-vector) perpendicular to the plane of incidence and "n" to 
polarization in the plane of incidence. Since skylight is 
also polarized, the equation above becomes two equations: 


L r " r l s l 


L « r s 
r ii n 


(id 


For small angles, i.e. directions within about 20° of nadir, 
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Sky luminance distributions 
in sun-zenith plane for 
several locations and solar 
zenith angles, e s . Data 
compiled by Austin [Ref. 3] . 
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the reflectivity is given by 


= r = ( \ .L __ a \ = 2 

1 \ n w + V 


1 % 


( 12 ) 


** 


the well-known Fresnel formula for dielectric reflection at 

normal incidence in which n = 1.341 is the refractive index 

w 

of seawater and n = 1.0 is that of air. 

a 


]' At other angles r is given by formulas that are straightforward 

*’ but complicated and are given in numerous physics and optics 

: r texts. They are of little if any practical value here because 

I. the sea surface is ruffled by wind which averages the reflectivity 

over a range of angles. Gordon [Ref. 4] has computed the average 

j 

I for various wind speeds with the results given in Fig. 18. 


I The wind-ruffled sea surface blurs the reflected image of the 

4 - 

sun into a glittering pattern called sun glint. The radiometer 

- 

} should be pointed away from the sun to avoid this glint as 

much as possible, Fig. 16. However, to reduce sky contamination 
J* to a minimum, the radiometer should look somewhat farther from 

..i 

the reflected sun angle, as Fig. 19 shows, about 75° when the 
sun is high in the sky to 90° when it is low. 


■i . 



i 


By coincidence, this same direction gives the maximum polarization 
of reflected skylight in the perpendicular direction due to both 
f • rs in Eg. 11: 


r jL 


> r ( and > S. 


ii 


The former fact is apparent from Fig. 17 and the latter results 
from the dipole radiation pattern of molecular scatter shown 
in Fig. 16. Unfortunately Fig. 18 is computed for unpolarized 
light and we do not have a similar one for the separate components. 
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When the sky is mostly clear, so that light has the directional 
qualities shown in Fig. 16, then the radiometer should be pointed 
to avoid reflected skylight and a polarization analyzer should be 
used to reject the perpendicular component. The radiometer's 
field of view should be centered about a direction 75 to 90° 
away from the reflected sun image in the plane containing that 
image and the nadir. This technique gives four advantages: 

Avoids sun glint 

Rejects sky contamination by looking at the darkest 
part 

Further rejects skylight by taking advantage of its 
polarization 

Takes further advantage of the polarization by 
Fresnel reflection. 

If the sky is clear (so that molecular scattering dominates 
over particle scattering) , then a polarization filter almost 
completely eliminates sky reflection at the longer wavelengths.* 
This is fortunate, because it normally causes most of the 
veiling luminance at the longer wavelengths; compare and 
L in Figs. 8, 9, and 10. 

w 

♦Some blue and violet remain, which have been depolarized by 

double scattering in the atmosphere, but this falls off as 
-8 -4 

X , a factor of a for each of two Rayleigh scatters. This 
is why the sky looks violet in some photos taken through 
Polaroid at right angles to the sun. 


> 

* 
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1 . 5 Transmission 

We shall use transmission factors T to account for loss of 
light in propagating through lossy media: 

T - atmosphere (satellite case) 

T - sea surface 
s 

Tp - polluted layer 

In every case, light that passes down through a layer must 
return through the same layer to be detected, and so we use 
capital T to represent round-trip propagation. However, this 
does not imply that the loss is the same each way. For example, 
.when the sun is rather far off the zenith, Fig. 20, and the 
radiometer is looking straight down, the transmittance through 
some layer going down is T = e a while going up it 

is considerably less, namely T = e aZa ; hence, 

T = TT = exp [— ot z a (1 + sec 0)] (13) 

In estimating T, light that is scattered through small angles 
should not be counted as lost. A cutoff angle is not easy to 
define, and will have to be determined by empirical matches 
to the remote sensing equations using data taken at sea. Fig. 21 
shows a summary of estimates of the unscattered transmittance of 
the atmosphere needed for the case of a satellite sensor. 


1.5.1 


Polluting Layer 


For measuring Vt^ a suitable absorption cell is shown in Fig. 22. 
The inside cylindrical wall is silvered for the same reason 
discussed in connection with Fig. 13. Similarly, the known 
incident irradiance overfills the stop aperture, and the 
detector is large enough to catch all light at the output. 


\ J radiometer 
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Fig. 21. Vertical beam transmittance of the atmospheric 
models (curves A, B, C, D) and as measured 
] by Guttman (curves a and b) . Applied Optics 7 

, (1968) p. 2377 
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Transmission varies approximately exponentially with distance z, 
H = e“ az 

but not exactly owing to multiple scattering in thick layers. 

The attenuation rate may be expressed 

a = a. + a C 

WO 

where a,, is the attenuation in the water that dilutes the 
pollutant, a Q is attenuation of undiluted pollutant, and 
C is its concentration. Combining these equations gives 

T = e~ awZ e ~ a °^ Cz ) (14) 

The factor Cz in parentheses is the area concentration of 
pollutant, so many grams per square meter, which says in effect 
that it makes no difference how the pollutant is distributed 
in depth. This is true so long as there is no chemical change 
that depends on concentration, as there is in the case of the 
DuPont acid dump. In the latter case, a Q is a function of C. 

It is prudent not to assume that Eq. 14 holds over any widv 
range of C or z until verified in the laboratory. Over the . 
intervals where it does hold, Eq. 14 provides an interpolation 
formula. 


1.5.2 


Surface Transmission 


Near vertical transmission through the sea surface is given by 

2 / „ \ 2 


( 4n n \ z / n \ 

(<) 


(15) 


The first factor is round-trip Fresnel reflection at the dielectric 

2 

interface ((1-r) , see Eq. 12). The second factor represents the 
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reduction of radiance rising out of the water due to angular 
spreading. Fig. 1. Using n fl = 1.0, n^ = 1.34 as before gives 

T s = (0.96) (0.56) = 53% (16) 

The fact that the first factor has little effect implies that 

surface roughness is not a big factor in determining T , as 

s 

shown in Fig. 18, except that patches of foam cause significant 
changes by scattering light upward and sideways. 
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2.0 


REMOTE SENSING EQUATIONS 


This part of the report derives equations for determining, 
pollutant concentration in two cases. The first is for a 
low altitude aircraft, and the second for high altitude, most 
likely a satellite. 

2.1 Low Altitude 

The significance of this case is that low altitude eliminates 
all atmosphere effects except reflected sky radiance, and we ‘ 
can deal with that by quite effective means. As discussed in 
Sec. 1.4, a good choice of aspect angle plus a polarization 
filter eliminates most of the interference from the sky when 
it is clear. At other times it is necessary to measure the 
spectral radiance of the sky directly with a sky radiometer 
so that the reflected radiance may be subtracted out of the 
sea data, Fig. 23. The best arrangement is to combine sea 
and sky radiometers into one instrument with the same detector 
and amplifier, so that any fluctuations in sensitivity affect 
both measurements equally. An optical switching or scanning 
arrangement causes the detector to look at the sky periodically 
between looks at the sea. (It need not look at the sky very 
often.) In any case, the pointing of the sky and sea radiometers 
must be at supplementary angles to the nadir as indicated by 
the reflection geometry in Fig. 22- 

When the sea has normal roughness, the reflected image of the 
sky is blurred. This causes no problem when the sky is uniform, 
either clear or overcast. However, when the sky has broken 
clouds, the optics in the sky radiometer must blur the direct 
image in a manner that simulates the blurred reflected image 
that the sea radiometer sees. The proper amount of blur may 
be determined by measuring the angular diameters of the eliptical 
sun glint pattern. 
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The radiance seen over clean seawater is given by 


L = HT R + rS (17) 

C 5 W 

As Fig. 23 shows, the first term represents light that penetrates 
the water, and emerges, and the second term light that reflects 
from the surface. The latter is L r , Eq. 11, either r,, S,, or 
r i|S|| + r _L s J_ depending on whether a polarization filter is 
used. The first term in Eq. 17. is the larger by a factor 
of about three through most of the visible spectrum, but the 
two cross over near the red end. Figs. 8 through 10. Recall 
that every quantity in Eq. 17 is a spectrum. Their colors 
are as follows : 



white 


R - blue to yreen (water color) 
w 

(blue when clear 
S ~ \ white when overcast 


Also recall that Eq. 17 implies three spectral measurements of 

light intensity. The radiance of the sea and sky, L and S 

c 

just discussed, and the irradiance H. 


Now consider the changes that occur when the water is polluted. 
An oil slick changes r considerably, but this discussion is 
restricted to the case of soluble pollutants, which have 
negligible effect on r and merely change the color of the. 
water. Fig. 24 shows the reflection and transmission factors 
that lead to the following equation: 

L = HT ( R + T R ) + rS (18) 

p s p p w 

The radiometer flying over a plume of pollution measures both 

L and L while also monitoring H and S. With these results, 
c p 
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Eqs. 17 and 18 are two equations to be solved for two unknowns. 

One unknown is the desired pollutant concentration implicit 

in its reflectance R and transmittance T . The other unknown 

P P 

can be one of the quantities R (water color) , T (surface 

w s 

transmission) , or rS (reflected sky) , that is not easy to 
measure or estimate accurately. There are different cases 
in which any one of these may be poorly known. However, we 
postpone that discussion to complete the mathematics of the 
remote sensing equations. 

When the color of the water R is unknown and eliminated 

w 

from Eqs. 17 and 18, the result is Eq. 19 in the table on 

the following page. When surface transmission T is unknown, 

the result is Eq. 20. This equation has the added advantage 

that irradiance H is eliminated and therefore does not have 

to be measured. When reflected sky rS is unknown, the result 

is Eq. 21. In general one must solve the appropriate remote 

sensing equation 19, 20, or 21, by an iterative algorithm: 

assume a pollutant concentration, test the corresponding 

R (A) and T (A) in the equation, and evaluate the fit. Assume 
P P 

another concentration and evaluate, and so on to find a best fit. 

A limiting case of interest occurs when the pollutant is a 

clear colored liquid, i.e. R * 0. In this case the remote 

P 

sensing equations (19 - 21) become Eqs. 22 and 23 in the table. 

Whenever > L c , this case cannot apply because the transmittance 
Tp has to be less than one. Eq. 22 is particularly nice because 
so many quantities drop out and therefore need not be measured. 

Another limiting case is T = 1, i.e. the pollutant does not 
appreciably shade the water beneath it.* This gives Eqs. 24 and 25. 

*Stri:tly speaking, if T =1, then R =0, but reflectances 

P P 

at sea are very small, so R can be .3% per steradian and T = 95% 

P P 

is close enough to 1. 
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TABLE OF REMOTE SENSING EQUATIONS 
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( 25 ) 





Equation 24 is a nice case because so few quantities are 
needed. 

Let us return to the postponed question , which quantity should 
be treated as the second unknown? Much experimental work is 
needed to establish the circumstances reliably, but some 
general guidelines are as follows: 

When the sea is fairly calm treat (water color) as unknown 
and use Eg. 19 . Sec. 1.3.2 discussed the fact that the color 
of seawater is not as obvious as it seems, and in fact measurements 
in situ are required to reliably determine 1^. Accordingly.. 

R^ should be treated as the unknown whenever the sea is calm 
enough to be free of whitecaps so that -T and r are fairly 

B 

well known. 

When the sea is choppy with foam and spray and the sky is 
uniform, assume T_ is unknown and use Eq. 20. One must work 
in a location where R^ has previously been established or else 
use only two colors, green and infrared, as discussed in 
Sec. 1.3.2. Foam and spray make it impossible to estimate 
either T or r. However, r may be measured indirectly leaving 
T as the unknown. For X ~ 700 nm, T is negligible, so Eq. 17 

S S as 

reduces to 

L ( A > 700) = rS ( A > 700) 
c 

The sky radiometer measured S directly, and so r(A > 700) = 

L /S. Moreover, the dependence of r on A is very weak and 
c 

known, so r is established throughout the spectrum leaving 

T as the unknown, 
s 
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When the sky and sea are too complex to find a reasonable 
est imate of blurred reflected sky radiance, treat rS as the 
unknown and use Eq. 21. When both T_ and rS are indeterminate, 

1 — - — — 11 1 " 1 * ' ' S 

rS should be treated as the unknown in the equations, because 

2 

T is never too far from the factor (n /n) fcr the spread 
s a v 

in solid angle. Other losses at the surface are not large 
(but not negligible either) . 


2.2 High-Altitude Sensor, Especially Satellite 

Let and denote the spectral radiance (clean and polluted) 
measured by a sensor at high altitude. These quantities are 
related to the low-altitude measurements discussed in Sec. 2.1 
as follows: 


L c m L c T a + T a 


L T + L* 
pa a 


(26) 


Here 7 • is the one-way upward transparency of the air path 

a 

and L* the radiance (blue) of the air along the path. In 
other words, the luminance we wish to observe is attenuated 
by T and veiled by L*. Unfortunately the path radiance L* 
is about seven times larger than the terms with the desired 
information. This means that L* must be subtracted off with 
considerable precision to avoid contaminating the ocean data 
with its errors. 


Recall that and in these equations are given byEqs. 
17 and 18: 


L c * HT s R w + rS (17) 

L p . HT S (R P * Vw> trS (18) 
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In Jec. 2.1 we processed these two equations in three different 
ways Appropriate to three different cases. However, in tha 
high- u‘ ! t it ude case, we always take the difference . 

- £ p - £ c - H Va (R p + Vw - V < 27 > 

This not only eliminates the largest term L* (and its uncertainty) , 
but also rS in which neither r nor S is known very well from 
high altitude. (Sea state and low-lying hydrometeors may be 
in doubt.; 

The factor H in Eq. 27 is irradiance at the sea surface, a 

quantity not normally available when the only measuring platform 

is at high altitude. Therefore, we replace it by the solar 

irradiance H measured at high altitude times some estimate 
s 

T a of downward atmospheric transmission: 

H = H s T r 

. Ai = H s (T a Ta) T s <R p + T p R w - R w ) 

It may seem that the product Ta T should be replaced by T 

3 l d 

to represent round-trip propagation as in Eq. 13 and other 
places. However, there is a subtle distinction between the 
two directions that needs to be emphasi ad. In the downward 
direction scattered light counts as long as it reaches the 
surface and thereby contributes to illumination. In the 
upward direction, this is not the case. As in Fig. 25, the 
polluted area is usually quite small compared with the altitude 
of the sensor. As a result the scattered light is diluted 
and completely lost in the scattered light from surrounding 
clean water. A distinction between polluted and clean areas 
is possible only in the unscattered light because the radiometer 
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images them and compares radiance on each side of the image 
boundary. Thus light scattered in upward propagation dx'ops 
out of the difference (Eq. 27) but unscattered light contributes, 
and we emphasize this with a chage of notation: 

D - T including scattered light (D stands for down 

3 3 . 

and diffuse) 

U - T excluding scattered light (U stands for up 
S3 

and unscattered) 

Finally we have the remote sensing equation: 

“ - H sWs l R p - d - V V (28 > 

The evaluation of this expression is not nearly as satisfactory 
as the corresponding equation for low altitudes (choice of 
Eq. 19, 20, or 21). The direct measurements of H (at the 
surface) and S (sky) in the low altitude case have been 
replaced by quantites that involve indirect evidence and even 
conjecture. Clouds near the line of sight affect D . Foam, 
spray, and possibly other low-lying hydrometeors affect T . 

One can estimate a correction based on whiteness of the reflection 
from clear areas. However, it cannot be very accurate because 
one does not know how much of the whiteness to attribute to 
clouds (reduce D ), surface foam and spray (reduces T ), or 
to reflected sky (clouds), which has dropped out of Eq. 28. 
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